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Metagenomic analysis was applied to bacterial communities associated with the shoots of two field-grown rice cultivars, 
Nipponbare and Kasalath. In both cultivars, shoot microbiomes were dominated by Alphaproteobacteria (51-52%), 
Actinobacteria (11-15%), Gammaproteobacteria (9-10%), and Betaproteobacteria (4-10%). Compared with other rice micro- 
biomes (root, rhizosphere, and phyllosphere) in public databases, the shoot microbiomes harbored abundant genes for CI 
compound metabolism and 1-aminocyclopropane-l-carboxylate catabolism, but fewer genes for indole-3 -acetic acid produc- 
tion and nitrogen fixation. Salicylate hydroxylase was detected in all microbiomes, except the rhizosphere. These genomic 
features facilitate understanding of plant-microbe interactions and biogeochemical metabolism in rice shoots. 
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MICROBES AND 
ENVIRONMENTS 



Rice plants are inhabited by complex microbial commu- 
nities that play important roles in plant growth and health 
(1, 18). Numerous ecological studies utilizing culture- 
dependent and -independent analyses have been conducted to 
investigate these rice-associated microbes (3, 14). These 
studies revealed the diversity as well as functions of rice- 
associated microbes. However, our understanding of micro- 
bial communities in rice is still limited with respect to certain 
groups of culturable microbes and their characteristics. 
Recent advances in high-throughput sequencers have enabled 
metagenomic studies to be conducted on various micro- 
biomes in different environments in order to investigate their 
phylogenetic and functional gene diversities. However, until 
recently, no metagenomic analysis had been applied to 
plant-associated endophytes due to the technical difficulties 
associated with preparing metagenomic microbial DNA that 
was not seriously contaminated with plant DNA. A technique 
to enrich bacterial cells from plant tissues was recently devel- 
oped (7) and was shown to be useful for metagenomic analy- 
ses of rice-root-associated microbiomes including bacterial 
endophytes (8). In the present study, we applied metagenomic 
analysis to microbiomes associated with the rice shoots 
(aerial parts) of two rice cultivars, Nipponbare and Kasalath, 
which were grown in a paddy field. 

Nipponbare (Oryza sativa L. japonica-type) and Kasalath (O. 
sativa L. indica-type) seeds were placed on two layers of filter paper 
in a petri dish (6-cm diameter) on 28 April 2008 and incubated at 
30°C. After 2 d, the germinated seeds were sown in a commercial 
soil mixture (Mitsui-Toatsu No. 3; Mitsui Toatsu Chemicals, Tokyo, 
Japan) in a 60-cm x 30-cm 800-cell tray (cell diameter, 1.5 cm; 
depth, 3 cm) and grown in a greenhouse under natural light condi- 
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tions for 4 weeks. The seedlings were transplanted into an experi- 
mental field at Tohoku University (Miyagi, Japan) on 29 May 2008. 
Rows were spaced 30 cm apart. Rice was cultivated under water- 
logged conditions with a water depth of 30 cm. Nine whole rice 
plants of each cultivar were carefully dug out by hand on 4 August 
2008. The roots were roughly washed by irrigation water to remove 
large soil particles, then immediately transported to the laboratory 
and washed with tap water to remove tightly adhering soil particles. 
The shoots were separated from the roots and stored at -80°C until 
they were used for subsequent analyses. The shoot samples were 
manually ground into a fine powder in liquid nitrogen using a mortar 
and pestle to release endophytic bacterial cells within the rice tissue. 

Shoot samples from the 9 plants were pooled and the meta- 
genomic DNA of each cultivar was extracted using a previously 
described cell enrichment method (7, 8). DNA samples were then 
subjected to metagenomic DNA sequencing. Shotgun sequencing 
was performed using a 454 GS FLX Titanium pyrosequencing 
system (Roche Diagnotics K.K., Tokyo, Japan). Pyrosequencing 
generated 1.2 million reads containing 479 million bases for 
Nipponbare and 0.8 million reads containing 342 million bases for 
Kasalath (Table SI). Sequence data have been deposited in DDBJ 
under project accession number DRA000376. Metagenomic reads 
were mapped to the rice genome sequence (IRGSP release, build 
5.0) using GS REFERENCE MAPPER (Roche Diagnotics K.K.) 
to eliminate contaminated reads from the rice genome (Table SI). 
Artificial replicate reads were removed using the 454 Replicate 
Filter (5) with default parameters (Table SI). As a result, 0.8 million 
reads containing 328 million bases remained for Nipponbare and 0.6 
million reads containing 240 million bases remained for Kasalath 
(Table SI), and these were used for later analyses. Phylogenetic 
analysis of metagenomic reads was conducted using MG-RAST 
(1 1) with the best hit classification mode of BLASTX. The maxi- 
mum E-value was set to 10~ 10 . GenBank was used for the database 
of BLASTX comparisons. Functional gene classifications of the 
metagenomic reads were conducted under the same conditions using 
the KEGG database. The relative abundances of functional genes 
were calculated by dividing the abundance by the total number of 
reads and multiplying the result by 1,000,000 (reads per million 
reads). 

Phylogenetic analysis revealed that the rice shoot micro- 
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biomes of Nipponbare and Kasalath were dominated by 
Alphaproteobacteria (51.1-52.4%), Actinobacteria (11.3- 
14.5%), Gammaproteobacteria (8.6-9.7%), and Betapro- 
teobacteria (4.5-10.1%) (Table 1). The abundance of 
Betaproteobacteria in Kasalath (10.1%) was more than 
twice that in Nipponbare (4.5%). At the genus level, 
Methylobacterium was the most abundant in both Nipponbare 
and Kasalath (Table 1). The abundance of Methylobacterium 
in Nipponbare (15.2%) was higher than that in Kasalath 
(9.3%). These results clearly indicated that the community 
structures of rice-shoot-associated bacteria differed between 
Nipponbare and Kasalath, as was previously shown using 
ribosomal intergenic spacer analysis (RISA) conducted in the 
same fields and with the same rice cultivars (15). The results 
of the present study further indicated that this difference was 
mainly caused by the high abundances of Methylobacterium 
in Nipponbare and Betaproteobacteria in Kasalath. The 
methanotrophs Methylosinus (0.2-0.3%) and Methylocystis 
(0.2%) were detected in both Nipponbare and Kasalath 
(Table 1). Rice paddies are known to be a major source for 
the emission of methane. Methane generated in rice paddy 
soil is transported through the rice aerenchyma and released 
from the rice culm into the atmosphere (13). Therefore, it is 
conceivable that methanotrophs colonize the routes of meth- 
ane transportation such as the aerenchyma. However, recent 
studies reported that some plants produce methane aero- 
bically from the pectin in the leaves (9); therefore, the extent 
to which methanotrophs contribute to the oxidation of meth- 
ane from the roots needs to be examined in future studies. 
Methane emitted from leaves may be the substrate for 
rice-shoot-associated methanotrophs (6). It is also notewor- 
thy that methanogenic archaea were detected in both 
Nipponbare (0.1%) and Kasalath (0.1%) (Table 1). Further 
studies are necessary to reveal the activities of methanogenic 
and methanotrophic microorganisms living in the rice phyto- 
sphere. 

Plant hormones are signaling molecules that are important 
for regulating the growth, development, and defense reac- 
tions of plants (2). Endophytic microorganisms frequently 
have the ability to produce or degrade plant hormones (2). 
We used MG-RAST to explore the rice shoot microbiomes 
for genes involved in the metabolism of plant hormones, and 
compared them with those in other metagenomes from rice 
roots (8), rhizosphere soil (10), and the phyllosphere (the 
total above-ground surfaces of rice plants) (10) deposited in 
public databases (Table S2). 

The gene for salicylate hydroxylase was found in the 
metagenomes of rice shoots, roots, and the phyllosphere, but 



Table 1. Phylogenetic compositions of bacterial communities 
associated with rice shoots (% abundance) in the cultivars 
Nipponbare and Kasalath 
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not in the metagenome of rhizosphere soil (Table 2). 
Salicylate hydroxylase mediates the conversion of salicylate 
into catechol. Salicylate is a phenolic phytohormone involved 
in pathogenic resistance (2). These results suggested that 
plant-associated bacteria may regulate the defense responses 
of the host plant by controlling the metabolism of salicylates. 
However, rice is known to have a high basal level of endo- 
genous salicylic acid (17), which suggests that rice-shoot- 
associated bacteria may also consume salicylate as a nutrient 
source. 

The abundance of the gene for 1-aminocyclopropane-l- 
carboxylate (ACC) deaminase was higher in shoot micro- 
biomes than in the other microbiomes (Table 2). This enzyme 
is responsible for cleaving ACC — the ethylene precursor in 
plants — into a-ketobutyrate and ammonia, thereby reducing 
ethylene levels in the plant, and the activity of ACC deami- 
nase produced by plant-associated bacteria can interrupt 
the biosynthesis of ethylene by plants (4). Ethylene is a 
plant hormone that can modulate various plant processes such 
as growth, development, and stress responses (20). Shoot- 
associated microorganisms may control ethylene levels to 
facilitate symbiosis as endophytes. Alternatively, ammonia 
produced by the hydrolysis of ACC may serve as a nitrogen 
source for shoot-associated microorganisms (4). 

Indole-3 -acetic acid (IAA) regulates plant growth and 
development. IAA-producing microorganisms have been 
isolated from rice tissues in previous studies (2). However, 
the gene for tryptophan 2-monooxygenase, which partici- 
pates in IAA synthesis via the indole-3 -acetamide pathway, 



Table 2. The relative abundances" of genes relevant to the metabolism of plant hormones 
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" Relative abundances are shown in reads per million reads. The KO number is based on the KEGG database (http://www.genome.jp/kegg/). See 

Table S2 in the supplementary materials for metagenome sequence accessions, rice cultivars, and sampling locations. 
b Root: DNA of the bacterial fraction (including endophytes) extracted from rice roots (8). 
0 Rhizosphere soil: DNA of rhizosphere soil around rice roots (10). 

d Phyllosphere: DNA of epiphytic bacteria released from the surface of rice shoots by sonication and shaking (10). 
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was not detected in the microbiomes of rice shoots (Table 2). 
In addition, the relative abundance of the indolepyruvate 
decarboxylase gene, a key enzyme in the indole-3 -pyruvic 
acid pathway of IAA biosynthesis, was lower in the shoot 
metagenomes than in other metagenomes (Table 2). These 
results suggested that IAA metabolism may be less important 
for plant-microbe interactions in rice shoots than in other 
parts of the rice plant. 

CI compound metabolism is prominent in Methylo- 
bacterium (16). The high abundances of Methylobacterium 
in shoot microbiomes (16) suggest that CI compounds such 
as methanol are major carbon and energy sources for 
shoot-associated bacteria. Methanol is first oxidized to 
formaldehyde by methanol dehydrogenase. The relative 
abundances of mxaFI and mdhl/2, which encode methanol 
dehydrogenase subunits, in shoot microbiomes were 5- to 
17-fold those in the other microbiomes (Table 3). Our data 
included endophytes within rice shoot tissues, whereas the 
reference data for the rice phyllosphere were derived from an 
epiphytic bacterial community on the shoot surface. Thus, the 
high abundance of methanol dehydrogenase genes implies 
that methanol is a crucial source of carbon and energy for the 
endophytic bacterial community in rice shoots. The reads for 
methanol dehydrogenase genes in rice shoot microbiomes 
were similar to those of Methylobacterium (47.8-51.9%), 
Rhodopseudomonas (21.7-27.5%), and Granulibacter (7.6— 
10.8%) (Table S3). These genera may be the major contribu- 
tors to methanol oxidation in rice shoots. 

Formaldehyde is a metabolite at a branching point in the 
dissimilation and assimilation pathways of C 1 compounds 
(21). Three main pathways have been identified for the 
dissimilation of formaldehyde: the tetrahydromethanopterin 



(H 4 MPT)-dependent oxidation pathway, glutathione-dependent 
oxidation pathway, and nicotinamide adenine dinucleotide 
(NAD)-dependent oxidation pathway (19, 21). Genes rele- 
vant to the H 4 MPT pathway were found in shoot microbiomes 
with abundances similar to those in other microbiomes (Table 
3). The abundances of genes for the glutathione-dependent 
oxidation pathway were more than twice as high in rice 
shoots and the phyllosphere than in roots and rhizosphere 
soil. The abundances of genes for the NAD-dependent oxida- 
tion pathway were higher in shoot and phyllosphere micro- 
biomes than in the others (Table 3). Formate resulting from 
the oxidation of formaldehyde is further oxidized to carbon 
dioxide, a reaction that is catalyzed by formate dehydro- 
genase (FDH). The major subunit of FDH was more than 
twice as high in rice shoot metagenomes than in roots and 
rhizosphere soil (Table 3), suggesting that CI compounds are 
important energy sources for bacteria in rice shoots and 
phyllosphere. Two main pathways have been identified for 
the assimilation of formaldehyde: the serine pathway and 
ribulose monophosphate (RuMP) pathway. The abundances 
of genes for the serine pathway were higher than those for the 
RuMP pathway in shoot metagenomes (Table 3), suggesting 
that the serine pathway is a major route for the assimilation of 
CI compounds. 

Nitrogen is also vital for microorganisms. The abundances 
of the genes for assimilatory nitrate reductase (nasA) and 
nitrite reductase (nirA) were higher in shoots and the 
phyllosphere than in roots and rhizosphere soil (Table 4), 
suggesting that those genes play important roles in the acqui- 
sition of nitrogen in shoot microbiomes. The abundance of 
narGHI, which encodes dissimilatory nitrate reductase, was 
lower in shoots and the phyllosphere than in roots and rhizo- 



Table 3. The relative abundances" of genes relevant to CI metabolism 
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NAD-dependent oxidation pathway genes 
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Formate dehydrogenase gene 
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a Relative abundances are shown in reads per million reads. KEGG database (http://www.genome.jp/kegg/). See Table S2 in the supplementary 
materials for metagenome sequence accessions, rice cultivars, and sampling locations. 
b Root: DNA of the bacterial fraction (including endophytes) extracted from rice roots (8). 
c Rhizosphere soil: DNA of rhizosphere soil around rice roots (10). 

d Phyllosphere: DNA of epiphytic bacteria released from the surface of rice shoots by sonication and shaking (10). 
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Table 4. The relative abundances" of genes relevant to nitrogen metabolism 
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1 Relative abundances are shown in reads per million reads. KEGG database (http://www.genome.jp/kegg/) 
b Root: DNA of the bacterial fraction (including endophytes) extracted from rice roots (8). 
c Rhizosphere soil: DNA of rhizosphere soil around rice roots (10). 

d Phyllosphere: DNA of epiphytic bacteria released from the surface of rice shoots by sonication and shaking (10). 



sphere soil, while the relative abundance of nirBD, which 
encodes dissimilatory nitrite reductase, was twice as high in 
the metagenomes of shoots and the phyllosphere than in roots 
and rhizosphere soil (Table 4). The abundances of nitrogen 
fixation genes (ni/HDK) in shoots and the phyllosphere were 
markedly lower than those in roots and the rhizosphere 
(Table 4). Nitrogen fixation may contribute to the fitness of 
microbes in the rice rhizosphere (10). The abundances of 
genes for nitronate monooxygenase (which oxidizes nitro- 
alkanes into the corresponding carbonyl compounds and 
nitrile) and formamidase (which hydroxylates formamide 
into formate and ammonia) were higher in shoots than in 
others (Table 4). These pathways may be important for nitro- 
gen acquisition in the endophytic bacterial community in rice 
shoots (12). 

The present study provides several notable insights into the 
metabolism of plant hormones, carbon, and nitrogen in 
rice-shoot-associated bacterial communities. Our results will 
assist in more deeply understanding plant-microbe interac- 
tions and will provide a useful set of data for metaproteome 
analyses (10). 
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